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Summary. — We present studies of double longitudinal spin asymmetries in semi- 
inclusive deep inelastic scattering using a new dedicated Monte Carlo generator, 
which includes quark intrinsic transverse momentum within the generalized parton 
model based on the fully differential cross section for the process. Additionally, 
we apply Bessel-weighting to the simulated events to extract transverse momentum 
dependent parton distribution functions and also discuss possible uncertainties due 
to kinematic correlation effects. 
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1. Fully differential SIDIS cross section 

The transverse momentum dependent (TMD) partonic distribution and fragmentation 
functions play a crucial role in measuring and interpreting information towards a true 3- 
dimcnsional imaging of the nuclcons. TMD PDFs can be accessed in several experiments, 
but the main source of information is semi-inclusive deep inelastic scattering (SIDIS) of 
polarized leptons off polarized nucleon. For SIDIS, the theoretical formalism is described 
in a series of papers [1, 2] using a tree level factorization [3] where the standard momentum 
convolution integral [4] relates the quark intrinsic transverse momentum in a nucleon to 
the transverse momentum of the produced hadron Put- 

In this work we present a model independent extraction of the ratio of polarized, gn,, 
and unpolarizcd, f lt TMD distributions using a Monte Carlo (MC) generator based on 
the fully differential cross section, in which we reconstruct the transverse momentum 
of the final hadron after MC integration over the quark intrinsic transverse momen- 
tum. In the MC generator we used the SIDIS cross section described in Ref. [1]. The 
Bessel-weighted asymmetry, providing access to the ratio of Fourier transforms of g±L 
and /i [5], has been extracted. The uncertainty of the extracted TMDs was estimated 
using unintegrated, transverse momentum non-factorized distribution and fragmentation 
functions. 
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The MC generator software employs the general-purpose, self-adapting MC event 
generator Foam [6] for drawing random points according to an arbitrary, user-defined 
distribution in the n-dimcnsional space. 

We consider the following SIDIS process 

(1) £{l) + N{P)^£(l') + h(P h )+X, 

were t -is the scattered lepton, N is the proton target and h-is the observed hadron 
(four- momenta notations given in parentheses). The virtual photon momentum is defined 
q = I — I' and its virtuality Q 2 = —q 2 . Following the Trento conventions [7], we use the 
virtual photon-nucleon center of mass system, where the virtual photon momentum q is 
along the z axis and the proton momentum P is in the opposite direction. The detected 
hadron h has momentum Ph. In the parton model the virtual photon scatters off an 
on-shell quark. The initial quark momentum k and scattered quark momentum k! have 
the same intrinsic transverse momentum component k± with respect to the z axis. The 
initial quark has a fraction x of the proton's light-cone momentum, while the produced 
hadron h has a light-cone momentum fraction z and transverse momenta component pj_ 
with respect to the scattered quark's momentum k' (see Ref. [I]). The fully differential 
cross section used in the MC generator is then given by [I] : 

= K{x,y)J{x,Q 2 ,k 1 _) x 



dxdydzd 2 p±d 2 'k 



(2) 



X] fl,q{ X i k ±) D hq{ Z 'P±) + Wl ~ e 2 gi L ^ q {x,k 1 _)D liq {z,p 1 _) 



where the summation runs over the quark flavors in the target. The kinematic factors 
K(x, y) and e, and the Jacobian J(x, Q 2 , fcjj arc defined in [I]. A represents the product 
of lepton and nucleon helicities (A = ±1). fi. q {x, k±) and g\L, q {x, k±) are the quark TMD 
distributions and D\ tq (z,p±) is the TMD fragmentation function for an unpolarizcd quark 

q- 

In many phenomenological studies of the semi-inclusive deep inelastic scattering, the 
transverse momentum dependence of distribution and fragmentation functions is fac- 
torized from the light-cone momentum dependence x and z, respectively. In our MC 
generator we use the Bessel-weighting method [5] to alleviate this approximations in 
the extraction of TMD distributions from SIDIS measurements. We discuss the follow- 
ing modified Gaussian (MG) expressions for the TMD distributions and fragmentation 
function, inspired by AdS/QCD [8, 9], in which x and k± (z and p±) are not factorized: 



(3) fi(x,k±) = /i(ar)-r- 2 , gi L (x,k ± ) = g 1L (x) 



x{l-x)(kl) fl ' — ' ^ — 'x{l-x){kl) gi 



(4) 



£>i(*,P±) = £>i(js) 



z(l-z)(p 2 ± ) 



TMD distributions non-factorized in x and k± are also suggested by the diquark spec- 
tator model [10] and the NJL model [11]. For a test extraction, the x and z dependence 
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Fig. 1. — (Color online) The effective (k^_(x)) in the modified Gaussian (MG) ansatz for the TMD distributions 
(see Eq. 3) as function of x for two different input parameters < k^_) obtained from MC events simulated within 
the range 0.50 < z < 0.52 and using 6 GcV (solid symbols) and 11 GeV (empty symbols) lepton beam energies. 



in Eqs. (3-4) we use the parametrizations f\{x) = (1 — x) 3 x 1 * 313 ) gn,(x) — fi(x)x - 7 
and D t (z) = 0.8(1 - zf , with widths: {k\) fl = 0.54GeV 2 [8], (k 2 ± ) gi = 0.8{k 2 ± ) fl 
and (p 2 ± ) = 0.14 GeV 2 . Fig. 1 illustrates the average effective transverse momentum 
(kj_) = (A; 2 _)/ 1 (l — (x))(x) as function of x from MC events for 6 GeV (solid symbols) 
and 11 GeV (empty symbols) incoming electron beam energies and produced hadrons 
within 0.5 < z < 0.52. Note that the value for (hj_) fa obtained from the MC events is 
always smaller than the implemented value due to energy and momentum conservation 
[12]. A similar non-flat dependence of the average quark transverse momenta versus x 
and z is also observed in NJL-jet model [11]. 



2. Bessel-weighted extraction of the double spin asymmetry All 

We present the extraction of the double spin asymmetry All, defined as the ratio 
of the difference and the sum of electroproduction cross sections for antiparallel, <r + , 
and parallel, a~ , configurations of lepton and nucleon spins, using the Bessel-weighting 
procedure described in [5] and applied in [13]. Within this approach, one can extract the 
Fourier transform of the double spin asymmetry, A J L ^ hT PhT \bT) , defined as 
(5) 

A Mb T P hT ) (b ) = v + (b T )~a-(b T ) = Qll (Pt) = /p^ E g gifo z 2 b 2 T )D\{z, b 2 T ) 
LL [T ' d+{b T ) + d-{b T ) auu{b T ) V E J2 q f?(x,z 2 b 2 T )Dl{z,b 2 T y 
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using measured double spin asymmetries as functions of Put [14], for fixed x, y, and z 
bins. Here bx is the Fourier conjugate of the Put- The Fourier transforms of the helicity 
dependent cross sections, ct ± (&t), can be extracted by integration (analytic models) or 
summation (for data and MC) over the hadronic transverse momentum, weighted by a 
Bessel function Jo, 

N ± 

(6) ~a ± {b T ) ~ S± = ^ MbTPhT,i) ■ 
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Fig. 2. — (Color online) Bcssol- weighted asymmetry vs bT with and without the correction, together with 
analytical and numerical comparison from the MC with PDFs and FFs from Eq. 3. See the text for more 
details. 



The Bessel-weighted asymmetry obtained from the simulated events is shown in Fig. 2 
as function of br with filled (red) circles, while the analytic expression f 1 , zbT ^ using 

fl(x,zb T ) ° 

(k'j_}g 1 and (kj_) f ± from the fits to k\ distributions from the same MC sample is depicted 
by the (blue) full line. For values &t < 6 GeV -1 , which corresponds to about 1 /to, 
the Bessel-weighted asymmetries could be extracted with an accuracy of 2.5%, although 
with a systematic shift. This clear systematic shift between the extracted and calculated 
asymmetries is due to the kinematic restrictions introduced by energy and momentum 
conservation^) as well as binning effects, which deform the Gaussian shapes of the k± 



( x ) In the light-cone coordinate system the transverse component is less or equal to the mo- 
mentum component along the light-cone vector. 
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and p± distributions. In experiments, there is always a cutoff at high PhT due to accep- 
tance and the small cross section, as well as a cutoff at small PhT where the azimuthal 
angles are not well defined due to the experimental resolution. These restrictions in P^t 
directly affect the extracted k± and p± distributions and yield to the mentioned distor- 
tion of Gaussian shapes which result in the systematic shift between the extracted and 
the calculated asymmetries. Obviously, this shift depends on experimentally introduced 
restrictions for the accessible PhT range. 

We discuss two approaches which take these conditions into account, one corrects 
the data (simulated MC events) the other applies limits on the integration range for the 
intrinsic transverse parton momenta when calculating the asymmetry. In order to correct 
the data, the contributions from the missing P h T ranges outside the accessible values is 
estimated by choosing a certain model for the parton transverse momentum dependence, 
e.g. a Gaussian distribution. This model dependent correction of data is shown in Fig. 2 
by the (blue) filled squares. The extracted asymmetry does now match the theoretical 
curve for values 6t < 6 GeV" 1 . Alternatively, limited numerical integration over intrinsic 
transverse momenta in the calculation of the asymmetry, where the integration limits 
correspond to the accessible experimental PhT range, yields a calculated asymmetry that 
describes correctly the experimental situation without introducing a model dependence. 
This is shown in Fig. 2 by the (black) open squares. 

In summary, we used a fully differential Monte Carlo event generator to study how 
the implemented quark transverse momentum distributions k± and p± are being changed 
due to the energy and momentum conservation low and kinematic correlations within 
the generalized parton model. We then use the obtained k± and pj_ distributions to 
apply the Bessel-weighting strategy for the extraction of TMD distributions. As example 
we calculate the Bcsscl- weighted double spin asymmetry A J ^ T PhT \bT) ■ Within the 
b T range of 6 T < 6 GeV" 1 , which corresponds to about 1 fm, the Bessel- weighted 
asymmetries could be extracted within 2.5% accuracy. 
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